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Abstract 

In this paper, we investigate the density-dependent incompressible nematic liquid crys¬ 
tal flows in n{n = 2 or 3) dimensional bounded domain. More precisely, we obtain the local 
existence and uniqueness of the solutions when the viscosity coefficient of fluid depends on 
density. Moreover, we establish blowup criterions for the regularity of the strong solutions 
in dimension two and three respectively. In particular, we build a blowup criterion just 
in terms of the gradient of density if the initial direction field satisfies some geometric 
configuration. For these results, the initial density needs not be strictly positive. 

Keywords: incompressible nematic liquid crystal flows; density-dependent; strong solu¬ 
tions; blowup criterion 


1 Introduction 

Nematic liquid crystals contain a large number of elongated, rod-like molecules and possess 
the same orientational order. The continuum theory of liquid crystals due to Ericksen and 
Leslie was developed around 1960’s [H, [ 2 I (see also j^). Since then, numerous researchers have 
obtained some important developments for liquid crystals not only in theory but also in the 
application. When the fluid containing nematic liquid crystal materials is at rest, we have the 
well-known Ossen-Frank theory for static nematic liquid crystals, see the pioneering work by 
Hardt-Lin-Kinderlehrer Q] on the analysis of energy minimal conhgurations of nematic liquid 
crystals. Generally speaking, the motion of fluid always takes place. The so-called Ericksen- 
Leslie system is a macroscopic descriptions of the time evolution of the materials under the 
influence of both the flow velocity held and the macroscopic description of the microscopic 
orientation conhguration of rod-like liquid crystals. In this paper, we investigate the motion 

^Email: gaojcl998@163.com, Tel:+8613450423982 
^Email: taoq060@126.com 
^Email: mcsyao@mail.sysu.edu.cn 




J.C.Gao, Q.Tao, Z.A.Yao 


of incompressible nematic liquid crystal flows, which are described by the following simplihed 
version of the Ericksen-Leslie equations: 

' pt + div(pM) = 0, 

_ (pu). + div(,„ 0 u) - div(2,(p)fl(„)) + VP = -Adiv(Vd 0 Vd). 

divM = 0, 

Jt + u-Vd = e{M+\Vd\^d), 

in n X (0,+cxd), where is a bounded domain with smooth boundary in M"'(n = 2 or 3). 
Here p, m, P and d denote the unknown density, velocity, pressure and macroscopic average of 
the nematic liquid crystal orientation respectively. D{u) = ^ is the deformation tensor, 

where Vm presents the gradient matric of u and its transpose. /i>0, A>0,6*>0 are 
viscosity of fluid, competition between kinetic and potential energy, and microscopic elastic 
relaxation time respectively. The viscosity coefficient p = p{p) is a general function of density 
and be assumed to satisfy: 

/iGC^[0, cxo) and p > 0 on [0, cxo). (1.2) 


Without loss of generality, both A and 6 are normalized to 1. The symbol Vd © Vd denotes 
the n X n matrix whose (i, j)—th entry is given by Vdj • Vdj, for i,j = 1,2, To complete 
the equation fll.ll) - fll.2p . we consider an initial boundary value problem for fll.ll) - fll.2p with the 
following initial and boundary conditions 


{p,u,d)\t=o = {po,Uo,do) 


u = 0, 


with 

dd 

du 


|do| = 1, divMo = 0 in fl; 
= 0 on 


(1.3) 

(1.4) 


where u is the unit outward normal vector to dQ. 

When the fluid is the homogeneous case, the systems fll.ll) - fll.4p are the simplihed model 
of nematic liquid crystals with constant density. When the term |Vdpd be replaced by the 
Ginzburg-Laudan type approximation term d, Lin hrst derived a simplihed Ericksen- 
Leslie equations modeling the liquid crystal hows in 1989. Later, Lin and Liu made 

some important analytic studies, such as the existence of weak/stron^olutions and the partial 
regularity of suitable solutions. Recently, Dai, Qing and Schonbek |8| studied the large time 
behavior of solutions and gave the decay rate in the whole space in with small initial data. 
Grasselli and Wu also considered the long-time behavior and obtained the estimates on the 
convergence rate for nematic liquid crystal hows with external force. They also showed the 
existence of global strong solutions provided that either the viscosity is large enough or the 
initial datum is closed to a given equilibrium. As for the case of |Vd|^(i, Huang and Wang 


10l | established a blow up criterion for the short time classical solutions in dimensions two and 


three. Recently, Li ll| proved the local well-posedness of mild solutions with initial data, 
in particular, the initial energy may be inhnite. 

When the huid is non-homogeneous case, we would like to point out that the system (II.R - 
fll.4p include two important equations, which have attracted large number of analysts’ interests: 


2 
























Strong Solution to the Density-dependent Incompressible Nematic Liquid Crystal Flows 


(i) When c? is a constant and /i is a function depending only on the density p, then system fll.lD - 
fll.dp are the well-known Navier-Stokes equations with density-dependent viscosity coefficient. 
First, Lions 12|] established the global existence for the weak solutions for the case of positive 
initial density. As for the uniqueness, Lions tolled us the fact that sufficiently smooth solutions 
are unique and any weak solution must be equal to a strong one if the latter exists. Later, Cho 
and Kim 13] proved the local existence of unique strong solutions for all initial data satisfying 
a natural compatibility condition for the case of initial density being not be strictly positive. 
He also built the following blowup criterion: 


sup (||Vp(f)||L. ||VM(t)||L2) = CX), 

0<t<T* 


(1.5) 


if T* is the maximal existence time of the local strong solutions and T* < oo. 

(ii)When p is a constant, the systems fll.ll) - fll.4p are a density-dependent incompressible hydro- 
dynamic flow of liquid crystals. When the term iVdpcf be replaced by the Ginzburg-Laudan 


type approximation term the glob al existence of weak solution is obtained in 14-16 

for each e > 0. Recently, Hu and Wu [l7| proved the decay of the velocity held for arbitrary 
large regular initial data with the initial density being away from vacuum in two dimensional 
domain with smooth boundary. As for the case of | Vdpd, Wen and Ding 1^ obtained the local 
existence and uniqueness of strong solutions to the Dirichlet problem in bounded domain with 
initial density being allowed to have vacuum. Since the strong solutions of a harmonic map 
can blow up in hnite time 19], one cannot expect to get a global strong solution with general 
initial data. Therefore, many researchers attempt to obtain global strong solutions under some 
additional assumptions. Wen and Ding [1^ also established the global existence and uniqueness 
of solutions for two dimensional case if the initial density is away from vacuum and the initial 
data is of small norm. Global existence of strong solutions with small initial data to three 
dimensional liquid crystal equations are obtained by Li and Wang in 20] for constant density 
case, Li and Wang in 2l] for nonconstant but positive density case, and Ding, Huang and Xia 


m 


22 


Recently, Li proved the global existence and uniqueness of strong solutions with initial 


data being of small norm for the dimension two and three in bounded domain in 23] and the 
initial direction held satisfying some geometric structure for the two dimensional whole space 


m 


24 


In this paper, we investigate the density-dependent incompressible nematic liquid crystal 
hows when the viscosity coefficient is a function of the density of the huid. More precisely, 
we establish local unique strong solutions to fll.ip - fll.4p . Then, we consider the possible break¬ 
down of regularity for the strong solutions. Firstly, We build up a blowup criterion for the 
three dimensional bounded domain with smooth boundary. Secondly, applying a logarithmic 
inequality, we improve the preceding blowup criterion by omitting the velocity in a two dimen¬ 
sional bounded domain. Lastly, if the initial direction satishes some geometric conhguration, 
we establish a blowup criterion just in terms of the gradient of the density in two dimensional 
space. 

Before stating our main result, we hrst explain the notations and conventions used through- 
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out this paper. We denote 



fdx. 


Let 

/:=/* + «• V/ 

represents the material derivative of /. For 1 < q < oo and integer k > 0, the standard Sobolev 
spaces are denoted by: 


Wo’" = {n e = 0 on ifo^ = Wq’^ 


For two nx n matrices M = {Mij),N = {Nij), we denote the scalar product between M and 
N by 

n 

M:N=J2 

*j=i 

Finally, we recall the definition on the weak space which is defined as follows: 


4 |/ G : ll/lliP = sup t\{xeQ: \ f{x) > t\}\p < cx)| . 
Now, we state our first result as follows. 


Theorem 1.1. Let Vt be a bounded smooth domain in M"'(n = 2, 3) and q G (n, oo) be a fixed 
constant. Suppose that the initial data {po,Uo,do) satisfies the regularity conditions 

0 < po £ Uq E Hq n H^, do E and |(io| = 1, 


and the compatibility condition 


— diy{2p{po)D{uo)) + VPq + div(V(io © Vdo) = \fpog and divuo = 0 in 12 (1.6) 

for {Po,g) E X . Then there exist a positive time Tq > 0 and a unique strong solution 
{p,u,d,P) of (ll.ll) - (ll.4l) such that 

0<pE C([0,To]; W^’"), pt E ^([O, To]; L''), 
u E C([0, To]; n H^) n L^O, To; W^’^), 

Ut E T^(0, To; Hi), y/pUt E T°°(0, Tq; T^), 

PeC'([0,To];Pi)nT2(0,To;W'0, 

dEC{[0,To];H^)nL^{0,To;H^), |d| = 1 in 

dt E C{[0,To];H^)nL\0,To;H^), du E L^{0,To;L^), 

for some r with n < r < min{q, and Qtq = 12 x [0, To]. 

After having the Theorem 11.11 at hand, we will build the following blowup criterion of 
possible breakdown of local strong solution for the initial boundary problem (ll.ip - fll.dh in three 
dimensional space. 
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Theorem 1.2. Suppose the dimension n = 3 and all the assumptions in Theorem li.il are 
satisfied. Let {p,u,d,P) be a strong solution of the initial boundary problem fll.ip - fll.4p andT* 
be the maximal time of existence. If 0 < T* < oo, then 

(ll Vp||l°o(0,T;L9) + II'w|Il'>i(0,T;L;i) + IIL“2(0,T;Lf?)) = (l-^) 

where ri and Si satisfy 

2 3 

-1-<1, 3 < Tj < oo, i = 1, 2. (1.8) 

Si Ti 

As a corollary of Theorem II .21 we will establish a blowup criterion for the density-dependent 
incompressible flow when the viscosity depends on the density in three dimensional domain. 
More precisely, if d is a constant vector, then we have the following corollary. 

Corollary 1.3. Suppose d be a unit constant vector and all the assumptions in Theorem \1.S\ 
are satisfied. Let {p,u,d,P) be a strong solution of the initial boundary problem fll.ip - fll.4p and 
T* be the maximal time of existence. If 0 < T* < oo, then 

(||Vp||l°°( 0,T;L-J) + llwlUyO.T;!,-)) = OO, (1.9) 

where r and s satisfy 

2 3 

—I— <1, 3 < r < oo. (1-10) 

S T 

Remark 1.4. The criterion foru in (11.9p is given by a Serrin type and is more general than 
the blowup criterion fll.Sp . 

Our next work is to improve the proceeding blowup criterion (11.71) by utilizing a logarithmic 
inequality for the two dimensional space, i.e.. 

Theorem 1.5. Suppose the dimension n = 2 and all the assumptions in Theorem li.il are 
satisfied. Let {p,u,d,P) be a strong solution of the initial boundary problem fll.ip - fll.4p andT* 
be the maximal time of existence. If 0 < T* < oo, then 

(ll^PlU“(0,T;Lq + ||Vd||L»(0,r;L-)) = OO, (1-11) 

where r and s satisfy 

2 2 

—I— <1, 2 < r < oo. (1-12) 

s r 

As a corollary of Theorem 11.51 we provide a blowup criterion for the density-dependent 
incompressible flow when the viscosity depends on the density in two dimensional domain, 
then we have the following corollary. 

Corollary 1.6. Suppose d be a unit constant vector and all the assumptions in Theorem \1.5\ 
are satisfied. Let {p,u,d,P) be a strong solution of the initial boundary problem fll.ip - fll.4p and 
T* be the maximal time of existence. If 0 < T* < oo, then 

^limJ|Vp||Loo(o,T;Lq = oo. (1.13) 
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Our next work concentrate on building blowup criterion the same as fll.lSp if the initial 
direction held satishes some special geometric conhguration. More precisely, we have 

Corollary 1.7. For any i{i = 1,2), suppose the i—th component of initial direction field doi 
satisfies the condition 


0 < dgi — *^0* < 1 or — 1 < doi < —dg^ < 0, 

where dg^ is defined in fIb.Sp . and all the assumptions in Theorem \1.5\ are satisfied. Let 
(p, u, d, P) be a strong solution of the initial boundary problem fll.ip - fll.4l) and T* be the maximal 
time of existence. If 0 < T* < oo, then 


lim ||Vp||l° 

2 ^_ yfP* 


“(0,T;L'J) 


= OO. 


(1,14) 


The rest of the paper is organized as follows: In Section 2, we present some useful lemmas 
which will play an important role in this paper; In Section 3, we prove the Theorem 11.11 by 
applying the method in 13|]; From Section 4 to Section 6, we discuss and verify blowup criterions 
of strong solution respectively. 


2 Preliminaries 

In this section, we give some useful lemmas which will be used frequently in this paper. The 
hrst lemma is the regularity estimates for the stationary Stokes equations, i.e.. 

Lemma 2.1. (See fl^J lAssume /i G C^[0, oo) and p G 0 < p < C. Let {u, P) G Hi x Lf 
be the unique weak solution to the boundary value problem 

—div{2p{p)D{u)) + VP = P, divM = 0 in 12; J Pdx = 0, 

where D{u) = YidiVT, Then we have the following regularity results: 

(1) //P G Lfi then (n, P) e x and 

W'^Wh^ + II-P|Ihi < c'II-^IIl 2 (i + ||Vp||l<j)'j-"'. ( 2 . 1 ) 

( 2 ) If F E for some r G (n, q) then (n, P) G x IF^’^ and 

_ qr 

+ ||-P||wr’’ ^ C'||P||Lr(l + IIVpIIl?). (2-2) 

(3) //P G H\ then (n, P) e x and 

||w||_H-3 + ||P||//2 < C||P||hi( 1 + llpllvr^'i)^ (2.3) 

for some N = N{n,q) > 0. The constant C depends also on \\d‘^p/dp^\\c. 

Next, we introduce a Holder inequality in Lorentz space. The Lorentz space and its norm 
are denoted, respectively, by and || • Wlp.i, where 1 < p < oo and 1 < g < oo. Now, we can 
state the following Holder inequality in Lorentz space 
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Lemma 2.2. (See \2al )Let 1 < pi,p 2 < oo with ^ ^ ^ ^ Qiy '?2 < C) 0 . Then for 

f E and g E it holds that 

11/ ■ S'IIlp,'? < C\\f\\LPim llghpi’ii with q = min{gi, ga}, 


where C is a positive constant depending only on pi,P 2 , Qi and q 2 - 


The following lemma has been proved in 2^, for the readers’ convenience, we give the proof 
in detail. 


Lemma 2.3. Assume g E and f E with r E {n, oo], then f ■ g E Lf. Furthermore, for 
any e > 0 and r E (n, oo], we have 

ll/•^?lli.<^ll^?ll?.l + C(e)||/||el|^?|li., (2.4) 

where C{e) is a positive constant depending only on e, n, r and the domain f2. 

Proof. Applying the Lemma 12.21 it is easy to get 


ll/'^IU^ = ||/■^||L 2,2 < C|]/|]l.|]^||^^, 2 , where r e (n,oo]. (2.5) 

Next, we will show that 

r — n n 

||^||^^.2 < C'|]^||^||^|]])i, where r e (n,oo]. (2.6) 

Indeed, if r = oo, then fl2.6p holds on obviously. If r G (n, oo), then is a real interpolation 

2’-! 2r2 

space of L’'i -2 and , where ri, r 2 and r satisfy n<ri<r<r 2 <oo and ^ = — + —, then 


< ciisir 


< c 


^1- 


ri 

L2 


1 . 

2 




1 

2 


r — n 


n 

r 


r/i’ 


where we have used the Sobolev inequality. Therefore, combining fl2.5p with (12. 6 p gives (12. 4 p 
directly if we exploit the Cauchy inequality. □ 


The last lemma introduced in this section will be the following logarithmic Sobolev inequal¬ 
ity which plays an important role in the proof of the Lemma 15.21 Omitting the proof for brief, 
one can read [^ . 

Lemma 2.4. Let Q be a bounded smooth domain in and f E L‘^{s,t]H^ fl for 

q E (2, cxo). Then there exists a constant C depending only on q such that 


2 


0 <c 


+ l|/lli2(s,t;Hl) lll(e + \\f\\L'^{s,t-,W^^1)) 


where C depends only on q and O, but independent of s,t. 


(2.7) 
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3 Proof of Theorem 11.11 

In this section, we will only give the existence proof for the Theorem 11.11 since the uniqueness 
of the solution is easy to obtain by a standard argument (c.f.[l^). In order to solve the initial 
boundary problem fll.ip - fll.4p . we will split the proof into three parts. In part one, we will 
establish the global strong solution for some linearized systems. In part two, we prove the 
solution of the linearized systems converges to the original initial problem fll.ip - fll.4p in a local 
time for positive initial density. In part three, we verify Theorem 11.11 for the case of general 
initial density. 

3.1. Global existence for the linearized equations 
We consider the following linearized systems: 

Pt + v - yp = 

put + pv ■ Vm — div(2/i(p)Il(M)) -|- VP = —div/, (3.1) 

divM = 0, 

in (0, cxo) X n, where 2D{u) = Vm+ p = p{p) satisfies fll.2p and u is a known divergence-free 
vector field. Then we will state the main result in this subsection. 

Proposition 3.1. Assume that the data {po,uo,f) satisfies the regularity conditions: 

0 < Po e Uo e n and / e L°°(0, T; H^) n L2(0, T; fi G L‘^{0, T; H^) 

for some q with n < q < oo, and the compatibility condition 

— div(2/i(po)P(Mo)) + VPo -|- div/o = y/pog and divwo = 0 in hi, (3.2) 

for some {Po,g) E x . If in addition, v satisfies the regularity conditions 

V e L°°(0, T; n p2) n ^^(o, T; IT^’"), Vt G ^^(O, T; H^) and divw = 0 in 12, 

for some r with n < r < min{g, Then there exists a unique strong solution {p,u,P) to 

the initial boundary value problem fl3.ip . fll.2p - fll.4p such that 

pG C([0,T];iy^’''), Vm, P G C([0,T];P^) nL2(0,T;IT^O, 

Pi gC'([ 0,T];P), VpWieL°°( 0 ,T;L 2 ), e L\0,T-H^,). 

In order to prove Proposition 13.11 we will take by three steps: 

(1) In addition to the assumptions in Proposition 13.11 if suppose 

p G C^[0, cxo), po G po> d for some 5 > 0, 

then we give the prove of Proposition 13.1( 

(2) To remove the additional condition (13.4p . we need to derive some uniform estimates inde¬ 
pendent of S, ||po||w 2.9 and \\d‘^p/dp‘^\\c; 
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(3) Having the results of the proceeding two steps at hand, it is a standard argument to give 
the proof of Proposition 13.11 

Now, let us to begin our hrst step. Indeed, taking the method in 


13|, it is easy to get the 


following results. For a brief, we only state the results and omit the proof. 


Lemma 3.2. In addition to the hypotheses of Proposition \3.1l we assume that the condition 
(13.41) are satisfied. Then there exists a unique strong solution (p, u, P) to the initial boundary 
value problem (13.ip . (I1.2p - (ll.4p such that 

p e ^(0, T; n L2(0, T; pt G L“(0, T; 

u e ^([0, T]; Hi n H^) n T; Ut G L°°(0, T; L^) n T; Hi), 

P G L°°(0, T; H^) n T; H'^), 

where n < r < min{g, 

Thanks to the previous Lemma [221 there exists a unique strong solution (p, u, P) satisfying 
the regularity (13.3p . To remove the additional hypotheses (13.4p . we will derive some uniform 
estimates independent of 5, ||po|| 1 ^ 2.9 and ||(9^p/(9p^||c- 


Lemma 3.3. Suppose {p,u,P) be the strong solution to the problem (13.ip . (I1.2p - (ll.4p . then we 
have 


sup (IIpIIwi.'j + WpiWli + ||m||h 2 + ll^ll^i + \\\fpUt\\L^) 

0<t<T 


+ 




2 

14/2.’- 


+ 11 ^^ 411 ^ 2 ) dt < C, 


(3.5) 


where C independent of 6, ||po||iv 2.9 and \\d‘^p/dp‘^\\c and n < r < min{g, 


Proof. Step 1: We deduce from (I3.ip i by applying the characteristic method that 


llp(^)IU'> = IIpoIIl” for 0 < t < T, 1 < s < 00 . (3.6) 

Taking the gradient operator to (I3.1l) i . multiplying by g|Vp|'^“^Vp and integrating by parts, 
we obtain 

||Vp(t)||L9 < ||Vpo||L9exp (c 


\v{s)\\w^,rds 


which implies 


due to (I3.ip i . It 


\\dtp{t)\\L<i < C||n||H2||Vpo||L9exp J \\v{s)\\w^,rds ) , 
is easy to observe from (II.2p and (13.61) that 

C-^<p<C and |Vp|<C'|Vp|, 


(3.7) 


which will be used repeatedly. 
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Step 2: Multiplying fll.ip o by Ut and integrating over (0,t) x 12, we have 




{\ft\\^u\ + p\v\\Vu\\ut\ + \p'\\v\\Vp\\D{u)\‘^) dxdr (3.8) 


p{p)\D{u)\'^{t)dx + I I p\ut\'^dxdT 

<C + J\f\\\/u\dx + 

4 

= C + '^Iu. 

i=l 

To estimate the terms Iii(l < < 4), we will make use of the Young and Sobolev inequalities. 

^11 ^ C'(s) [ \f\‘^dx + e [\Vu\‘^dx, 


h2<2 (II/,Hi. + ||Vn||i.)dr, 

Jo 

h3<C{e) [ [ p\v\‘^\'Vdxdr+ e 



p\ut\‘^dxdT 


< C{e)\\p\\L^\\v\ 


Loo 




I'Vul'^dxdr + e / / p\ut\'^dxdT, 


h 4 <C / |1v|U6||Vp||l‘j||Vm|P 12, dr 
Jo 


rt 12q —n(q+6) 

< c / ||Vt,|U.||Vp|U,||v«||y •• ||v«||„. 


n(q+6) 

dr 


<C||Vp|U,||Vn|U2 / liVnIli.dr + e / ||Vn||^idr. 


On the other hand, we get 


p{p)\D{u)pdx > C / \D{u)pdx = 


2 C 


\Vu\‘^dx, 


(3.9) 


due to fl3.7p and fl3.1l) o. Substituting Ji/l < i < 4) and 03.91) into 03.Sp yields 


1 

40 


iVnl^da; + 



p\ut\‘^dxdT 


<C l+\\f\\l.+ \\ft\\UT)+C WVuWUr + e \\Vu\\Ut. 


(3.10) 


In order to deal with the term || VM||i.i(ir, we will applying the regularity estimate for the 
stationary Stokes equations in Lemma [2. II More precisely. 


where 


\ u \\ h 2 + IIPII^I < 0||F|U2(1 + llVpiuo^ < 
11^11^2 = II - put - pv -Vu- div/||i 2 

< 0 ||y^M,||i 2 + 0||Vm||l 2 + ||div/||i 2 . 
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Then we have the following regularity estimate 

\\u\\h^ + \\P\\m < C(||div/||i2 + \\^/put\\L^ + ||Vm||l2). 

Substituting fl3.1ip into (13.101) and choosing £ small enough, we obtain 


(3.11) 


\'Vu\‘^dx + 



p\ut\^<C + C \\Vu\\l,dT, 


which, if we exploit the Gronwall inequality, implies 


sup llVwllis + / W^UtWl^dr <C. 
0 <t<T Jo 


(3.12) 


Step 3:Differentiating fll.ip „ with respect to t, multiplying by Ut and integrating over fl, we 
have 


[ p\utfdx + ^ [ iVutfdx 


2 dt 


C 


<C j (|n||Vp||VM||VMi|+p|n||M4||VMt| + |nnVp||VM||Mi| 

5 

+p\vt\\Vu\\u,\ + |/,||V'lI,|)dl = Njii- 


(3.13) 


2=1 


To estimate the terms / 2 i(l < i < 5), we will apply (13.lip , the Gagliardo-Nirenberg and Holder 
inequalities repeatedly: 

hi < ||w||L->||Vp||L'i||VM||^^||VMt||L 2 

<G(e)||n||i.||Vp||L||V«||^3+£||V«,||i^ 

< G(e)(l + ||div /|||2 + \\^/pUt\\l 2 ) + e\\Vut\\l 2 , 

h2<C{e)\\p\\L<=^\\v\\\^ j p\ut\^dx + e J\Vut\^dx 

<C{e)\\^ut\\l2+e\\Vu,\\l,, 
hs < ||w|lioc||Vp||L,||VM||^^||MjL6 

Sg+nq —3n 3q+3n —nq 

<Giin|ii.iivpiu.iiv«ii^r^ iiv«ii;;7^ \\vu,h2 

< C{e){l + lldiv/lllz + \\y/put\\L2) + e\\'Vut\\l2, 

hi < ||p||l,oo||nJi3||VM||L2||Mt||i6 

< C\\p\\L^\\VVt\\LA\'^u\\L2\\VUt\\L2 

< C{e)\\'Vvt\\‘i2 + +e\\'Vut\\‘i2, 
hb < C'(^)ll/t|li 2 + ^||VMt||^ 2 . 

For any hxed r G (0, t), substituting / 2 i(l < * < 5) into fl3.13p and integrating over (r, t) C [0, T] 
yield 


2 /pi“.r*+^ 



{Vutfdxds 

rt 


< 


p{T)\ut{T)\'^dx + C (||div/||^2 + ||/t|li2 + ||Vnt||i2)cis+ p\utfdxds 
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Thanks to the compatibility condition fl3.2p . letting r —?■ O’*" and applying the Gronwall in- 
eqnality, it arrives at 

snp \\^/pUt{t)\\l 2 + [ \\Vut\\l 2 dt < C. (3.14) 

0<t<T Jo 

Step 4: High order estimates. Indeed, combing fl3.14p with fl3.lip - fl3.12p yields 

-l- ||T’||jyi < C. (3.15) 

Applying the stationary Stokes regularity estimate, i.e. fl2.2p . we get 

_ qr 

||'u||y[/2,r -|- ll^ll VKi’’’ — C||f|U.(l + ||Vp||„)’^<C'||F|U.., 

where 

ll-^IU"- = II - put - pv-Vu- div/llir 

< llplU°°ll“tlU’' + llplk°°ll'^lk°°l|VM||L'- + ||div/||ir 

< C'(IIpIU°°I|VmJl2 + IIpIU-II^IIl°°I|Vm||hi + ||div/||Hi) 

< + ||VMt ||£2 -|- ||div/||jyi). 

Hence we obtain the following regularity estimate 

11^11 14/2,r -f- ll^ll C{\ -\- ||V'Ut ||/,2 -|- ||div/||j^i). (3.16) 

Therefore, we complete the proof of lemma. □ 

After having the Lemmas 13.2113.31 at hand, we turn to prove the Proposition 13.11 We only 
sketch the proof here since it is a standard argument (c.f.[l3|). Let (po,uo) be an initial data 
satisfying the hypotheses of Proposition 13.11 For each S G (0,1), choose p^ G and G 
G^[0, oo) such that 

0 < (5 < Po < Po + 1) Po Po in and —)■ /i in cxo), 

as (5 —0, and denote {uq, Pq ) G Pq x a solution to the problem 

—div{p^{pl)D{uQ)) + WPq + div/o = and divug = 0 in fl. 

Then, applying the Lemma the corresponding solution {p^,u^,P^) satishes the estimate 
sup f||p'^||wi'9 + IIPiIU'J + ll'n'^||_H '2 -|- ||P'^||jyi -|- II llna'j 

0<t<T V '' 

+ [ (||'n‘^|1^2,r -|- llP'^lIwi-’’ A ll^'nj||^2) dt < C, 

Jo 

where C independent of 5, Upolliv^.? and ||5^p/5p^||c and n < r < min{g, We choose a 

subsequence of solutions {p^,u^) which converge to a limit (p, u) in a weak sense. Therefore, it 
is a strong solution to the linearized problem satisfying the regularity estimates in Lemma 13.31 
Thus, we complete the proof of Proposition 13.11 
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Strong Solution to the Density-dependent Incompressible Nematic Liquid Crystal Flows 


3.2. Local existence for the original problem 

In this section, we will prove a local existence result on strong solution with positive initial 
density to the original problem fll.ip - fll.4p at first. Furthermore, we derive some uniform bounds 
which are independent of the lower bounds of the initial density. Then, these uniform bounds 
will be used to prove the existence of strong solution with nonnegative initial density in the 
last part of this section. 

Proposition 3.4. Assume that the data {pOiUo^do) satisfies the regularity conditions 

PoefF'’", uo^HlnH\ doeH\ 
for some q with n < q < oo and the compatibility condition 


— diY{2 p{pq)D{uo))+ VPo + d\v{Vdo QVdo) = ^/p^g and divwo = 0 in hi, (3.17) 

for{Po,g)eH^ X L^. Assume further that po > S inQ for some constant 5 > 0. Then there exist 
a time To G (0,T) and a unique strong solution {p,u,P,d) to the nonlinear problem (ll.ip - (ll.4p 
such that 

peC'([0,To];lF^’''), AeC'([0,To];L''), 

u e C([0, To]; H] n 77^) n T2(0, To; IT^’^), 
ut e T2(0 ,To;77o'), G T°°(0, To; T^), 

FGC'([0,To];77i)nT2(0,To;lF'0, 

dGC([0,To];773)nT2(0,To;77"), |d| = 1 in 
dt G 7;([0,To];77i) nT2(0,To;772), du G L\0,To-,L^), 
for some r with n < r < min{g, 

To prove the Proposition 13.41 we first construct approximate solutions, as follows: 

(1) first define = 0 and dP = do] 

(2) assuming that and dP~^ was defined for 7 > 1, let (p^, dP, P^) be the unique solution 
to the following initial boundary value problem 

' pk + . vp^ = 0, 

p^Ui^ + - div(2p(p'=)71(u^)) + =-div(Vd^ 0 Vd^), 

^ T fc n (3.19) 

divM = 0, 

^ d’l - Ad’^ = - (u^-^ ■ V)d^-\ 

with the initial and boundary conditions 


(p^, d^) = (po. Mo, do) X G fl, 

. . ddK 


(U 


’ du 


■) = (0, 0) on 572, 


where u is the unit outward normal vector to 572. 


(3.20) 

(3.21) 
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3.2.1. Uniform bounds 


Thanks to the Proposition 13.II to fl3.191) ^ — fl3.19L and existence and uniqueness of the theory 
of parabolic equation to fl3.19p ^fc.f. |27l|). it is easy to get the existence of a global strong solution 
with the regularity fl3.18p to the linearized problem fl3.19p - fl3.2ip . 

From now on, we derive uniform bounds on the approximate solutions and then prove that 
the approximate solutions converge to a strong solution of the original nonlinear problem. Let 
P > 1 be a hxed large integer, and define a function as 


^xit) = max sup (l 

l<fc<Xo<s<t ^ 


l|V«‘(i>)||i. + ||V<i‘(s)||„. + ||Vp‘(i,)||„). 


Observe then that 


5<p^ <C, C-^ <11^ <C, iVp'^l < C\Wp^ 


(3.22) 


Then we will estimate each term of ^k{1) in terms of some integral of ‘Fir(t). 
Lemma 3.5. There exists a positive constant N = N{n,q) such that 

\\^u\t)\\l,+ f \\J^^u\s)\\l,ds<C + C f<!>K{s)^ds 
Jo ^ Jo 


(3.23) 


for all k, 1 < k < K. 

Proof. Multiplying fl3.19p „ by integrating by parts and making use of fl3.20|) ^. we have 


= - j p'^iu'^-^ ■ W)u^ ■ utdx - J ■ Vp'^)\D{u'^)fdx + j Vd'^Q : Vwfda;, 

On account of the identity 

J Wd^QWd’^ : Wu'ldx = ^J Vd^QVd^ : Vu^dx- J Wd^QWd^ : Vm^ +V^QV^ : Vu^dx. 

We integrate over (0, t) and apply fl3.9p and fId.lhP g to deduce that 
^ f iVu'^l^dx- j\Vd^\^\Vu^\dx + ]^ [ [p'^lu’^l^dxdT 


'0 


<C + C 



+ iVd^llVd^llVw''!) dxdr ( 3 . 24 ) 


C + hi- 


i=l 
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Strong Solution to the Density-dependent Incompressible Nematic Liquid Crystal Flows 


Applying the Gagliardo-Nirenberg, Holder and Young inequalities repeatedly, it arrives at 

hi<C [ <6- / \\Vu'^-^\\l4Vu'^{^\\Vu%,dT 

Jo Jo 

<e f \\^/f^u'l\\l 2 dT + C{e) [ 


>0 


>0 


h 2 <C \\Vp1\l4^u^ \\l4^u^\\ 12 . dr 

Jo 

rt 12 q — n{q+ 6 ) 

<c / \\Vp^\\l4^u’^~"\\l2\\Vu^\\^, 

Jo 

<e f \\s/pu'l\\l.jT + C{e) f <S>1?dT, 


II V« 


(.+ 6 ) 
fc|| 6<j 7 

1^1 ar 


/33< ||Vd'=|U6||V<|U2||VM"|U3dr YG j ||Vd'=||Hi||V<|U2||VMi5"l|VMl||idr 

<s [ \\y/^u’l\\l 2 dT + C{s) [ 

Jo Jo 

for some Ni = Ni{n, q) > 0(i = 2, 3,4), where we have used the following regularity estimate 
\\u^\\m + \\P’"\\m < C{1 + \\^/f^u^t\\L^)^K for some Ni = Ni{n,q) > 0. (3.25) 

Substituting I^i^i = 1,2,3) into fl3.24p and choosing £ small enough yield 


1 f i^^ki2i^..ki^^ , 1 f p^\u\\^dxdT <C+ C I (3.26) 


— / \Vu^\^dx- / \Vd'^\^\Vu^\dx+- 



2C 

for some = N^{n, q) > 0. In order to control the term — J | on left hand side 

of (I3.26p . taking V operator to (I3.19p then we have 


- VAd^ = V - (u^-^ ■ V)d^-^] . 


(3.27) 


Multiplying by 4| Vd^pVd^, integrating (by parts)over f2 and exploiting the boundary condition 
ddS 


du 


an 


= 0, we get 


^ f \Vd'^\^dx + A j\Vd^\^\M'^\^dx 


<C j + |VdY|Vd^"^|=^ ^3 33 ) 

5 

+ |VdY|Vd^"^||V^d''-^| + iVdYlV^d'^l^) dx = 5 ^/ 4 *, 
where we have used the fact 


i=l 


- / VAdY4|Vd''|"Vd'^dxc= / 4|Vd''r|Ad'=rda; + 8 / djd^d^didid^did^did^dx. 


15 












J.C.Gao, Q.Tao, Z.A.Yao 


Applying the Holder and Gagliardo-Nirenberg inequalities, we have 

/43 < ||Vd^||ie||Vd^-iie < C\\Vd%4Vd^-YH^ < 

/44 < ||Vd1|ia||Vd^-ii6||VV-ii3 < WVdYmW^d'^YlHAYd'^YlH^ < C^k- 

h, < livdliiociiv^d^iii. < CllVdYmlYdY^ < c^K- 


Substituting lY = 15 2,3,4,5) into fl3.28p and integrating over (0,t), it arrives at 
[ iVd’^Yx + A [ [\Vd^\^\Ad’^YxdT <C+ C [ ^%{s)ds. 


Choosing a constant C* sufficiently large such that 

- ivw'^iivdY + aivdY > + ^|VdY, 

Z\_/ 40 Z 


then fl3.26p + fl3.29p x C* yields 

j\Vu'^Yx + j pYt?dxdT <C+ C <h^«dr, 
for some = N^n, q) > 0. Hence, we complete the proof of lemma. 


(3.29) 


□ 


Next, we estimate the term \\^^/p^u^\\l 2 and ||VMi||i 2 to guarantee the higher regularity. 
Lemma 3.6. There exists a positive constant N = N{n,q) such that 


\\Vf^u’l{t)\\l 2 + / llVnYx^ds <Cexp 
Jo 


C / <h^(s)ds 

Jo 


(3.30) 


for any k, 1 < k < K. 


Proof. Differentiating fl3.19p „ with respect to t, multiplying by and using fl3.19p ^. then we 
get 

^ ^ ^ p^Y\^dx + ^ f iVu’ll^dx 


2 dt 
< C 


C 



k\ k-hi^.k 
LL 




7 

+|c»,p‘l|Vt.‘’||V«f|)<ix = 5^/,.. 


(3.31) 


i=l 
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To estimate the term / 5 i(l < i < 7), we make use of the Holder, Gagliardo-Nirenberg and 
Young inequalities. 


<c||p'=|ihiiv«'=-‘iu.|iv^«fliyiv^«flli,liv«fiu. 

< c(E)\\i/^unh<i>i+ 

fe < l|plL-ll«‘-‘IU<l|Vt<'-|li.||v«‘|U.||«f|U. 
<C(£)||V«‘-‘ry||«*f„=+£||V«f||i. 

< C(t)(1 + \\F?ut\\l,)¥*^’^' +£||V«f|ry, 

/53 < l|p‘IU»l|v«*-‘||i.||v=«'=|U=||«f|U. 

<C{E)\\Vu>‘- 4 Ua'‘\\l ,+2 + E\\Vu>;\\l, 

< c(t)(i + |iV/«flli.)4Y"‘ +E|iv«flli^, 

/54<l|p‘iu»||v«*-‘|ii.||v«‘||„.||v«‘iu= 

<C(e)||Vu'=-‘||i,||«‘|||„+2 + £||Vuf||i. 

< ^(e)!! + +£||V«f||i., 

< iip'=ii!-ii«r'iu<iiv«‘iitfiivT«fiu» 

<c||v«f-'||y||v«'=iu.||y7«f|||.||yt*«f||j. 

< C(e,p)||Vp‘uflli>4>t+s||V«f|li.+ p||V«f-‘|||„ 
he < C(>;) /1V<i‘H+ £||Vt.fIII. 

<C(T)||V<i‘|||„||V4|||.+T||V«||||. 

<c(T)||v<i*|||,.||v4iii.+£||v«f|||., 

where we have used fl3.25l) . In order to control the term ||Vd^||L 2 , applying the estimate 
to the 03.271) . we obtain 

\\"^d^t\W < l|V (Ad^ + ■ Vd^-^) \\l2 

< C + llVd^-i^allV'd^-'lUa + ||Vd^-iie 

+ l|Vd^-iL^||Vn^-^|U2 + ||n^-i^e||VV-i^3) 

< 

Substituting 03.32p into Jse to deduce that 

h,<C{e)^], + e\\Vv!l\\l,. 
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Since the term J 57 is somewhat complicated, we will deal with it as follows: If n = 2, then 


h7< I \u’^-^\\Vp^\\Vu^\\Vu^\dx 


< IIVpIU^ I|v<||l2 

< C(e)(l + II ||i.)4"/“‘ + s|lV« 


fc ||2 

t Wl'^i 


If n = 3, then 


hi < ||Vp"|U,||M'=-ii6||VMi i.\\Vu\\\l 2 

£j Q d 

2q—6 g+6 

< ||Vp‘||„||Vii‘-‘|u.||Vti‘||j' ||Vti‘||i>i||Vtif|U« 

<C{e)1.^‘+e\\VuX^, 

where we have used the regularity estimate 

Wu^Ww^.r + WP^Ww^.r < C{1 + \\Vu'1\\l2)^^'' for some N7 = Nhn, q) > 0. 
Substituting / 5 j(l < i < 7) into fl3.3ip and choosing e small enough, we get 

1 


1 d 

2 dt 


+ — / ivufp* < (1 + liy/tiflii.)*? + >)llvuf-‘ 


- 1||2 


2C 


lL2- 


Fixing r in (0,T) and integrating over (r, t) C (0,T), we have 

llh/“J(i)llii + ^ llVafllii* 

< C||y/«f(r)|||, +cj\l + Iiy/«f||i.)tj"<ii> + lj‘ WVupYi^ds. 

From the recursive relation of ||Vmj ||l 2 , we get 

livufiii.*<cI -U) (|ivF«fWiii.+ ^‘(l + \\^/?A\\h)K’ds 


< 2C 


. 2 = 1 




\\\\ffhu'l{T)\\l2 + I {l + \\^/fhu’l\\l2)<^K^ds 


(3.33) 


Hence, we have the following estimate 

llVF«f(i)lli. + r liv«f|ii,* < 2 C { II vFuf(T)|ii, + /’(! + llvF^liO't?*') ■ 


Thanks to the compatibility condition (I3.17h . making use of the Gronwall inequality, it arrives 
at ^ ^ 

\c f S 2 ’(s)*l , 


ll\/F<4'(<)Ui= + / llViifllisds < Cexp 
Jo 

which completes the proof of the lemma. 


□ 
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As a corollary of Lemma 13.61 we can obtain the following estimate immediately. 
Lemma 3.7. There exists a positive constant N = N{n,q) such that 

||Vp^(t)||L'j < Cexp 

for any k, 1 < k < K. 


Cexp(C / $^(s)ds) 

Jo 


(3.34) 


Now we turn to the estimate the second term ||V(i^| 1/^-2 in ^k- Indeed, we should obtain 
the following estimate first. 


Lemma 3.8. There exists a positive constant N = N{n,q) such that 


f*t pt 

jk\\2 ^ yO I yO / , 


||Vd'^(t)||^i + / \\Vdt\\i2 ds<C + C / 

Jo Jo 

for any k, 1 < k < K. 

Proof. Step 1: Multiplying fl3.19p ^ by and integrating (by parts) over hi, we get 
1 d 


(3.35) 


2 dt 


\Vd^\^dx+ / \Ad^\^dx = / ■ V)d'=-^ - ■ Ad’^dx. 


Applying the Holder and Gagliardo-Nirenberg inequalities, it arrives at 

\jt /+ 

< J (|M^-^||Vd^"^||VV| + 

Integrating the proceeding inequality over (0,t), we have 

[ \Vd’^\^dx+ f [ lAd'^l^dxdr < C + C f 4>^dr. (3.36) 


Step 2:Multiplying fl3.27p by Vd^ and integrating (by parts) over hi yield 

- 4 - f |Ad¥dx+ f \Vd'f\^dx 
2dt J J 

< C{e) [ (|Vn'=-^nVd'=-^p + \u^-^\^\V^d^-^\^ + |Vd^-^|® + iVd^-^HV^d^-^p) dx 


+ e I \Vd^\^dx, 


19 











J.C.Gao, Q.Tao, Z.A.Yao 


where we have used the Cauchy inequality. Choosing e small enough, integrating over (0,t) 
and applying the Holder, Sobolev inequality, it arrives at 


\Ad^{t)\^dx + 



<C + C 



\Vd;\^dxds 


vv-y + dxds 


<c + c (II + II V«‘-‘||=.||VV-‘||(,. + ||Vrf*-‘||»,. 

Jo 

+l|vrf‘-‘|IPl|W-‘f„)ds 


<C + C I 


Hence, if applying the elliptic regularity, we obtain 

J\V^d'^{t)\^dx + j \Vd^^\^dxds<C+ C <h^ds, 

which, together with fl3.36l) . complete the proof of this lemma. 

Now, we can derive the high order estimate for d^. 

Lemma 3.9. There exists a positive constant N = N{n,q) such that 


□ 


||V=^d'=(t)||i.+ / \\V^dnhds<Cexp 


C / 


(3.37) 


for all k,l < k < K. 

Proof. Taking V operator to fl3.19p multiplying by VAd^ and integrating over fi, we have 


~ f \VAd’^\^dx+ j\Ad!l\^dx 


= j^J ■ Vd^-^ - \Vd^-^\‘^d’^-^) ■ VAd^dx 


(3.38) 


A 

dt 


[V(n^-^ • Vd^-^) - V(lVd^-^pd^-i)] ■ VAd^dx. 


Now we need to estimate the second term of right hand side of fl3.38l) as follows: 

^ ^ [V(n^-^ ■ Vd'^-^) - VdVd^-ipd^-^)] • VAd^dx 


dt 


< C / (iVurdlVd'^-d + iVn'^-dlVdf-d + lurdlV^d'^-d + |M''“d|VX 


TlVd^^-dlV^d^-d + |Vd^-d|VX^-d + |Vd''-d|V"d'=-d|d; 

8 


fc-l|IV7^fc-l| I L,fc-l||v72^fc-l| _|_ |^fc-l||^2^fc-l| 

k-1 1 I v 72 jfc-11 I jfc-11 


(3.39) 


+ |Vd^-ddVd^-d) iVAd'^lda; = ^hi. 


i=l 
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To estimate each term Iei{l < i < 8), applying Holder, Gagliardo-Nirenberg and Young in¬ 
equalities, we obtain 

hi < < C\\Vd’^-^\\H4^ut%4^^d%2 

<c<^% + c\\\/ut4h, 

h2 < \\Vu'^-4L4^dt4L4^"d%2 < c\\vu’^-4H4^"dt4L4^''d4L2 

< G(5)(l + + 6\\V^dt4h, 

hs < \\utX4^"d^~%4^''d%2<C\\Vut4L4'^"d’^~"M 

< + c\\\/u1-41d 

hA < \\u^- 4 L 4 \^^d\- 4 L 2 \\v^d 4 p 2 < c\\u^- 4 H 4 ^^d 4 L 2 \\v^(ii- 4 L 2 

< G(5)(l + + 5|| 

/65 < \\Vd’l-4L4V^d’^-4Ls\\V^d%2 < G||VX"-lLH|V'd"-irfi||V3d"|U2 

<C{6)^^j, + 6\\V^dt4lD 

he < \\vd>^-4L4\^"dt4L4^"d%2 < c\\vd>^-4H4^''d\\L2\\v^dt4L^ 

<C{6)^% + 6\\V^d’l-4l., 
hi < \\^d^-4L4\v^d^-4L3\\dt4L4^"d4L2 

< c\\wd>^-4H4^"d'^~"\\H4dt"\\H4^''d4L2 < 
les < ||Vd^-iie||VdriUa||V=^d^|U2 < 

<C{6)<^% + 6\\V^dt4h, 


where we have used fl3.25p . fl3.32p and fl3.33p . Substituting hi4 ^ i ^ 8) into fl3.39p . integrating 
fl3.38p over (0,t) and applying fl3.30p . then we have 


- / \V^d^\^dx + 





< 


j (|VM''-^||Vd'=-^| + + {Vd’^-^f) \VAd^\dx 


+ C exp 


C / 


^/7* + Gexp G [ +6 [ 

I Jo } Jo 


+ d / \\VX-^\\hds 
Jo 

+ 6 n\v^d'i-4i.ds 


(3.40) 


By Holder, Gagliardo-Nirenberg and Young inequalities, we obtain 

4—rt n 

Ill < \\^d^-4L4\vu^-4L4^^d^\\L^ < c\\vd’^-4j \\vd^-4^^4Vu’^-4L2\\v^d4L2 

4 —n n 

<c\\\/d^-4H4^u^-4L2\\v^d>^\\L2 + c\\vd>^-4^\\\/^d^-4l4\/u’^-4L4^^d>^\\L^ 


< e\\W^d4l,+v\\V^d^-4h+C{e,v) ( + \\Vd^-4l4^n^-Xl 


73 Jfc-l||2 


fc-l||2 


,fc-l||2 


7fc-l||2 


,k—l II 4 —r 
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2jk-l\ 


73 jfcl 


6 —n n 


<e\\W^d’^\\l+v\\V^d^-rL^+C{e,v) ( + || 


73 j/c-l||2 


12 

.k —1 II 6—n 


Its < ||Vrf"-iL6||V^rf"-ii3||V'rf1U2 < C\\Vd^-^HA\^^d'^-^J \\V^d 

12 —n n 

< C\\Vd’^-^\\l4V^d’^\\L2+C\\Vd’^-^^t \\^''d’^-^\\l,\\V^d%2 


72^A:-l||^^||V3^fc||^, 


2(12-n) 
-ik — 111 6—n 


< sWW^dYL^ + vW^^d^-Ti^ + C{s,v) + ||Vci'=-i^? 

/74 < \\vd^-m4v^d%2 < \\vd^-4%4v^d%2 
<e\\V^d4l.+Cie)\\\/d’^-4%,. 

Substituting J 7 j(l < i < 4) into fl3.40p . applying fl3.23l) and fl3.35l) and choosing e small enough, 
we obtain 


^ [\V^d^\^dx + 


<v\\^^d’^-4l, + s 



\V^d’l\^dxds 



\V^dr^\^dxds + C{6,r])exp 


C / 


Choosing 6 and rj suitably small, we have 


Iiv^d1li2+ / \\V^d4hds<i-{\\V-^d^-4h+ I \\V^dr\\hds]+Cexp 


3 jfc-l||2 


C / 


By recursive relation, it arrives at 


11 ^111.+ / ||VVf|| = .&<C j;— exp 


. 2 = 1 


C / <h^“ds 


which completes the proof. 

Now, from Lemmas 13.5113.91 we conclude that 


□ 


$x(t) < C exp 



Cexp ( C / 4>^ds 


for some N = N{n, q) > 0. Thanks to this integral inequality, we can easily show that there 
exists a time Tq G (0,T)(See Lemma 6 [281]), depending only on the parameters of C such that 

sup $x(t) < C. 

0 <i<To 

Therefore, using the previous Lemmas 13.5113.91 and other three estimates, we can derive the 
following uniform bounds: 

sup (iblk... + IbflU. + + IMIh. + livdfllid 


0<t<To 


rTo 


11^2 + ||dit||^2) dt < C, 


for all k > 1. 
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3.2.2. Convergence 

We next show that the whole sequence of approximate solution converges to a 

solution to the original problem fll.ip - fll.4p in a strong sense. To prove this, let us dehne: 

_fc+i _ pk+i _ -k+i ^ ^k+i _ ^k jk+i ^ ^k+i _ ^k^ 


Step 1; It follows from the linearized momentum equation fId.lbD o that 


pk+i^k+i ^ ^k+i^k . y-fe+i _ aiv[2p{p^^^)D{u^^^)] + V(P^+^ - P^) 




- p^+^u^ ■ - p^u’^ ■ + div [2(p(p^+^) - p{p^))D{u^)] 


(3.41) 


- div(Vd^+^ O O V(F+^). 

Hence multiplying fl3.41l) by and integrating (by parts) over fl, we get 


id //+‘|fi‘+‘|"<ij: + i- 


2 dt 


C 


\Vu 




dx 


= J {-d+'d - P+'k* ■ - p‘=u‘‘ ■ Vu* + div [2(^(p‘+‘) - p{p'‘))D(u'‘)] ,3 42 ) 

4 

-div(Vd^+^ O + V/ O Vd^+^)} • u’^+^dx = ^ hi. 


2=1 


To estimate hi{^ h i h 4:), using Holder, Sobolev and Young inequalities, we obtain 

hi < \\p^^^\\L4Ut - ■ VM^||L3||#+iL6 < C\\p^^^\\L2\\u^t - || l3 || ^2 

< C{e)\\p’^+4l4u1 - • Vu4ls + e||V#+ii 2 , 

h2 < \\p4U\\'/?A\l4^A\lA\u^^4l<i < \\p4U\\^/?u4l2\\wu4l4Wu^+4l2 

<C'(e)||/||i»|iy^s‘||i4|V«‘||i. + £||Vfi‘+'||i., 

Im<C j |y+‘||Vti‘||Vfi‘+‘|<i2'< CW j |y+y|V!i*p<i2: + e 

< c(£)iiv«*f„,,,iiy+‘iii.+£iivii‘+‘iii., 

h^ < 


|Vh 




dx 


I (|Vd^+^| + \Vd^\)\Vd^+ 4 Vu'^+^\dx 


<C{e) / (|Vd^+^|^ +iVd'^niVd'^+^l^dx + e / |Vh 


fc|2v 


7 fc+l| 2 ^ 


7-fc+l|2 


dx 


< C(e)(||Vd^+i^2 + ||Vdi^2)||Vd^+ii2 + +e |V#+^pdx. 


Substituting /^^(l < i < 4) into fl3.42p and choosing e small enough, it arrives at 

< c(||«f - u" ■ + ||v«‘|||,,,)|iy+‘|ii= + c||p'=iu»||v«‘|||,||v/7s'=lli. 

+ C4vd^+4l, + \\vd4lh\\vd’^^4l2. 
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Step 2: Since we have ■ V= 0, then multiplying by and 

integrating over f2 yield 

~ ^ - /«‘ • vp‘ ■ y+Vi 

<cj |p‘+y|V«‘|<ii + ||Vp‘||„||p‘+'||i.||fi‘||^^ (3.44) 

< C||v«‘iu«||p‘+'||'i. + ||Vp‘||„||p‘+‘|u.||vs‘IU. 

< C(i)(||Vp‘||„,.„ + ||Vp‘|||,)||p‘+'||i, + i||Vfi‘||i„ 

where we have used Young and Sobolev inequalities. 

Step 3: since we have 

- A(?+^ = • (Vd^ + Vd^-^) • + iVd'^-^pd*^ - (# • V)d^ - {u^-^ ■ V)d^ (3.45) 


Multiplying fl3.45p by —Ad^"*"^ and integrating (by parts) over f2, we obtain 



(Vd'^ + Vd^"^)d + iVd'^-^l^d^ - - (M^"^'^)d^] • (-Ad^+^)da: ^3 


2=1 


9^ 


To estimate / 9 i(l < i < 4), using Holder, Young and Sobolev inequalities, we have 


/91 < J |Vd^|(|Vd^| + |Vd^-i|)|d^||Ad^+Va^ 

<C{e) j |V<?p|V(i*| + |V(i*-‘|)^(ix + £||A<F+'||i, 

< C(£)(||V<i‘||i„ + ||V<i‘-'||i„)||V<i‘||i, +£||Aj‘+'||i., 

Iin<C{e) j\Vif-Y\f\^dx + e\\Af+'\\l, 

<C(£)||Vd‘-'||l.||<?||i. + £||Ad‘+'||i. 

< C(e)||Vrf*-‘||},.||VJ‘||i. +£||Aj‘+'||i.. 

/gg — — y djU^did^djd^'^^ + u^didjd^djdd'^^dx 

<C j (|VfiyV<i''||Vj‘+'| + |fi‘||VV||Vt?+‘|)<ix 

<C||V<i‘|U»||Vfi*||ig||Vj‘+'||ig+C||fi‘|U.||V"d‘||i,||V<i‘+'||ig 

<C(i)||Vd'=|||„||Vj‘+'|||,+CW||W||i.||V<F+'||i,+i||Vfi‘||i, 

< C(i)||V<i‘f„>||VJ‘+'||ig 

lM<C(e) j |M‘’-y|V(Fp(ix + £||A(i‘’+"|||g 
<C(£)||«‘-'|||„||V<F||i.<ix + £||Aj‘+'||ig. 
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Substituting Igi{l < i < 4) into fl3.46p and choosing e small enough, we get 

+ C{5)\\Wd^\\l4Vd^^^\\l2 + 5\\Vu'^\\l,. 


Denoting 




llL 2 

k(+\\\2 


+ Iivd 1 li., 




F^{t) = IK - n" • Vn'ilis + llVn'IlK^ + ||Vp"||i, + + ||Vd"+1^2, 


,k\\2 


k\\2 


jfc ||2 


jfc+l||2 


and 


jfc-l||4 


llVd'^+^l 


m- 


G^{t) = llpIlKllVn^lia + + ||Vd ii^a 

Then fiTTdIl + flT4T|l + fldlTIl yields 

K^^(t)+ [ it’^+\s)ds< [ [Cip\s) + 6it\s)] ds + [ C'K(s)K^^(s)ds, 

Jo Jo Jo 

which implies, by virtue of Gronwall inequality, that 

+ f 'ijj'^~^^{s)ds < C exp{Cst) f [C(p^{s) + ds. 

Jo Jo 

Choosing <5 > 0 and Ti so small that 8C(Ti + 5) < 1 and exp(C 5 Ti) < 2, then 




Kw + / i>''{s)ds 


for all 1 < k < K. By iteration, we have 


+ / tp'"^^{s)ds < 


^ 1 
fc+i/ ' ' - ■‘- 


qfc 


(p^{t)+ / 'ijj^{s)ds 


k>0. 


Together with the Poincare inequality and elliptic estimate, we get 

+ II '\/p*^h^’'"^|lL°o(0,T;I,2) + ||d^''’^||Loo(0,T;Hl) 


Up' 




|L°°( 0 ,r;L 2 ) 


Therefore, we get 


+ \\'u^~^^\\l^{o,t-,h^) + |M^’^^||L 2 (o,r;rf 2 ) < —. 

OO 

||P*''^^||l°°(0,T;L2) < OO, 

fe=l 

OO 

y~l ll'h^'''^||L 2 ( 0 ,r;_H'l) < C) 0 , 

k=l 

OO 

y^(IM^''"^IU°°(0,T;Hl) + ||d‘^’'"^||L2(0,T;rf2)) < OO, 


k=l 


(3.47) 


25 
















J.C.Gao, Q.Tao, Z.A.Yao 


which means 


as fc —)■ oo. 


i=2 

oo 

^ 

i=2 


1=2 


in L°°(0,T;L2), 
in L\0,T;H^), 
in L°°{0,T;H^)nL'^{0,T;H^), 


3.2.3. Conclusion 

Now it is a simple matter to check that (p, m, d) is a weak solution to the original problem 
fll.ip - fll. 4 p with positive initial density. Then, by virtue of the lower semi-continuity of norms, 
we deduce from the uniform bound that (p, m, P, d) satishes the following regularity estimate: 

sup (IIpIIwI'-j + WpiWli + \\^/pUt\\L^ + \\u\\h^ + \\p\\m + \\d\\m + llVdtHia) 

0<t<To 

+ [ (11^11^2,r + ||.P||^l,r + ||V^d||^2 + ||V^dt||^2 + 11 11 ^2 ) dt < C. 

Jo 

3.3. Proof of Theorem M.l\ 

Let (po,Mo,do) satishes the assumptions in Theorem 11.11 For each 5 > 0, let Po = Po + <5, 
satisfying 

Po ^ Po in as (5 —)■ O’*", 

and dg = do- Suppose {u^, P^) G Hi x is a solution to the problem 

—div(2p(po)T)(Mo)) + VPq^ -|- div(Vdo 0 Vdg) = \[^o9 uud divug = 0 in kl. 

Then by the regularity estimate, it is easy to get 

Mg G idg n and Mg —)■ uq in as <5 ^ O'*". 

Then, by Proposition l3.4l there exist a time Tg G (0, T) and a unique strong solution (p"^, P^, d^) 
in [0, Tg] xkl to the problem with the initial data replaced by (pg, Wg, dg). Note that (p^, P^, d^) 

will satishes the regularity estimate where C independent of the parameter 5. Hence, let 5 —)■ 0’*', 
it is a simple matter to see (p, w, P, d) is a strong solution and have the following regularity 
estimate 

sup (||p||wi.9 + llPtlU? -|- ||A/pMt||i2 -|- ||m||h2 -|- ||P||//i + ||d||//3 -|- ||Vdt||^2) 

0<t<To 

+ [ (11^11^2,r -f- ||P||^l,r + ||V^d||^2 -f" ||V^dt||^2 “f- 11 djt 11 ^2 ) dt < C. 

Jo 

Therefore, we complete the proof of Theorem 11.11 
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4 Proof of Theorem 11.2 


In this section, we will give the proof of Theorem 11.11 by contradiction. More precisely, let 
0 < T* < oo be the maximum time for the existence of strong solution (p, u, P, d) to fll.ip - fll.4p . 
Suppose that fll.7p were false, that is 

^0 = (^11 Vp||Loo(0,r;L'?) + \\M\l‘‘1(0,T;L^) + II^*^11 (0,T;LJ? )) < (^-l) 

Under the condition iu, one will extend existence time of the strong solution to fll.ip - fll.4p 
beyond T*, which contradicts with the dehnition of maximum existence time. 

Lemma 4.1. For any 0 < T <T*, it holds that 


sup 

Q<t<T 


L^ + \\^pu\\l 2 + \\Vd\\l,)+ j j (\Vd\^+\\Vd\^d + M\^^dxds<C, (4.2) 


where 1 < m < oo and in what follows, C denotes generic constants depending only on 
VL,Mq,T* and the initial data. 

Proof. Multiplying fll.ip ^ by < m < oo) and integrating the resulting equation over 

f2, then it is easy to deduce that 

||p(^)IU-= IIpoIIl- (l<m<oo). 

Multiplying fll.ip „ by u and integrating (by parts) over 12, it is easy to deduce 


(4.3) 


j p\u\^dx + — j \Wu\^dx < — j {u ■W)d ■ ^d dx. 

Multiplying (ll.ip ^ by Ad + | and integrating (by parts) over 12, we obtain 

|Vdpdx + J ||Vdpd +Adj^dx = J {u ■ \/)d ■ Addx + ^(dt + m ■ Vd)|Vdpd dx. (4.4) 

By virtue of |d| = 1, we have the fact 

(dt+ M-Vd)|Vd|2 d = 0. 

Substituting fl4.5l) into fl4.4p . it arrives at 

\jt I / ||V<i|^(i + A£ip£ia:= j (u ■ V)d ■ Addx, 


(4.5) 


which, together with fl4.3p . gives 


/(p|w|^ + |Vd|^)dx + ^ f |VM|^dx+ f I |Vd|^d + Ad|^ dx < 0. 


(4.6) 


Integrating fl4.6p over (0, t) yields 


2 / (Pl“l^ + |Vd|^)dx + 


which completes the proof. 



^|Vm|^ + I |Vd|^d + Ad|^^ dxdr J (po|moP + |Vdo|^) dx. 


□ 
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Now, we give the estimate for ||Vm||l 2 and ||V^d||L 2 . 

Lemma 4.2. Under the condition fl4.ip . it holds that for 0 < T < T*, 

sup (iiv«f„ + iivd|it. + iiv^dny 


0<t<T 


+ 


{\y/pu\^ + + |V^d|^ + \Vdt\^) dxdr < C. 


(4.7) 


Proof. Step 1: Multiplying fll.ip ^ by Ut and integrating (by parts) over fl, we have 


d 


d 


— / jj.{p)\D{u)fdx — — / Vd 0 Vd : Vudx + / p\ufdx 


dt 


dt 


= — j pit ■ {u ■ Vn) + p'{u ■ Vp)\D{u)\‘^ + Vdt 0 Vd : Vu + Vd 0 Vd* : Vudx. 
For any 0 < s < t < T, integrating over (s, t) and applying the Cauchy inequality yield 


(4.8) 


1 


I Vnl^da; — / Vd 0 Vd : Vudx + 



p|-u|^dzdr 


<C / |VMp(s)dx+ / \\/d\^\\/u\{s)dx + C{e,6) (plnplVn^ + |n||Vp||VM 





+ 1 Vd|^|V m|^) dzdr + e j j p|-u|^dzdr + d j j |Vdt|^da;dr. 
Choosing e = |, we obtain 


1 


\Vu\‘^dx — / VdOiVd : Vudx + 



p\u\‘^dxdT 


<C j (|Vm| 2 + |Vdn(s)dx + C(d) j j {p\u\^\Vu\^ + |n||Vp||Vn 
+ |VdnVMp)dxdr + d [ [\Vdt\^dxdT. 


(4.9) 


Estimate for the term J |n| | Vp| | Vnpdxdr. Indeed, by Cauchy inequality, Holeder inequality, 
and Sobolev inequality, we obtain 



|m| I Vp| I Vul^dxdr < C{e) 



<C{e) 

<C{e) 



|MnVMfdxdr + £ / |Vp| 


m|^| Vupdxdr + £ 

iv«r" 



iVpHVupdxdr 


Il§ 


9 dr 



|m|^|V ul^dxdr + £ / ||VM||^idr. 


In order to control the term ||VM||^idr, by virtue of the Lemma [2.11 we have 

\\u\\m + \\P\\m < C\\F\\l2{1 + ||Vp||lO^ 

< C|| — put — pu ■ Vu — div(Vd 0 Vd)||i 2 


(4.10) 


<C(||Vph|U2 + |||Vd||V2d|l 


L2>- 
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Hence we get the estimate 

J j \u\\Wp\\Wu\^dxdr 

J \u\^\VufdxdT + e J^{\\,/pu\\l 2 + \\\Vd\\V^d\\\l,)dT. 
Substituting fl4.1ip into fl4.9p and choosing e = it arrives at 


(4.11) 


\Vu\‘^dx — J Vd © Vd : Vudx + 



p\u\‘^dxdT 


2CJ - ./ • . ..-. 2.,, 

<C j {\Vu\^ + |Vd|^)(s)dx + C{d) j j {\u\^\Vu\^ + iVd^VMl" 

+ \Vd\‘^\V^d\‘^) dxdr + 6 J \Vdt\‘^dxdT. 

Step 2: Taking V operator to (inii4. then we have 

Vdt - VAd = -V(m ■ Vd) + VdVdpd). 

Multiplying fl4.13p by 4|VddVd and integrating (by parts) over 11, we obtain 

^ j |Vddda: + 4 J |VdnV^ddda: + 2 J\V{\Vd\^)\^dx 

= 2 [ I Vd|2 < VdVdp), z/ > da + 4 [ VdVd^d) : iVd^Vd dx 
Jon J 

3 


(4.12) 


(4.13) 


(4.14) 


4 / V(ti. Vd) : I V<i| W<i dx = ^ IIu, 

'' i=l 


where u is the unite outward normal vector to dfl. To estimate IIu = 2 | Vdp < V(|Vdd), z/ > 

da. Indeed, applying the Sobolev embedding inequality ly^’^(ll) L^{dfl), it is easy to get 


//,i<4 / |Vd|=>|VMd<T<C|||Vd|=>|VM||„„,. 

J dfl 

<C J dVdldV^dl + iVdldV^dp + |VdnV^d|)dx 

<C{p) [{\Vd\^\V^d\^+ \Vd\‘^+\Vd\^)dx + p [\V^d\^dx. 


(4.15) 


To estimate II 12 = 4 J V(|Vdpd) ; |VdpVd dx. Indeed, since |d| = 1, we have d ■ Vd = 0. 
Then, we get 

V(|Vd^d) : iVd^Vd = (V(|Vdd)d + IVd^Vd) : iVd^Vd = |Vd|®. 

Hence, it arrives at 

II 12 = 4 / V(|Vddd) : iVddVd dx = A [\Vd\^dx. (4.16) 
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By the Cauchy inequality, we have 

II 13 = -4 J V(m ■ Vd) : |Vrf| Vrf dx 

< [{\Vu\\Vd\^+\u\\Vd\^\V^d\)dx 


(4.17) 


<C J {\Vd\^\Vu\^ + iV^r + \u\^\V^d\^)dx. 

Substituting fl4.15l) - fl4.17l) into fl4.14p . choosing £ small enough, and integrating over (s,t), it 
arrives at 



J \\/d\^dx + A 

< j \VdWs)dx + Ciri) 

+ |Mp|V^dp) dxdr + r] 


\Vd\^\V^d\^dxdT + 2 




iVdVdHpdxdr 

(|Vdd + |Vdr + + iVdHVnp 

t /* 

V^d\^dxdT. 


(4.18) 



Step 3:Multiplying 04.131) by VAd, integrating (by parts) over hi and applying Young in¬ 
equality, we obtain 

1 d 


2 dt 


\Ad\^dx+ / \VAd\^dx 


< j (|VM||Vd| + InllV^dl + |Vd|^ + |Vd||V^d|)|VAd|da: 

<C{e) j{\Vu\^\Vd\^+ \u\^\V^d\^+ \Vd\^+ \Vd\^\V^d\^)dx + e j\VAd\^dx 
Choosing £ = | and integrating over (s,t), we get 


I Adl^dx + 



I VAd|^dzdr 



(4.19) 


< 2 / \Ad\^{s)dx + C / (|Vd|dV'dr + \^d\^\Vu\^ + |n|d+ \Vd\^)dxdT. 


Step 4:Multiplying (14.131) by Vd*, integrating (by parts) over 12 and applying Young in¬ 
equality, we obtain 

~ J lAdpdxY J |Vdiddx = y[V(|Vddd) - V(n-Vd)] ■ Vdtdx 

<C{e) J{\Vu\^\Vd\^+ \u\^\V^d\^+ \Vd\^+ \Vd\^\V^d\^)dx + e J\Vdt\^dx. 
Choosing £ = | and integrating over (s,t), we get 


lAdpdx + 



iVdiddxdr 



(4.20) 


< 2 / |Add(s)dz + C / (|Vd|dV"dd + \Vu\^\Vd\^ + \u\^\V^d\^ + |Vdr)da:dr, 
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Then, choosing 6 small enough and some constant C** suitably large such that 

- VdO Vd : Vm + + |Vd|^), 

ZLn 4(_/ 

then fl4.12p + fl4.18p x C** + (14.191) + fl4.20p and choosing r] and 6 small enough yield 


2m2 


J (|VMp + + \Vd\^){t)dx + J J {p\u\'^ + + |V(|Vd|")| 

+ \V^d\^ + \Vdt\^)dxdT 

<C + C [ (IVwp + + |Vd|^)(s)da; + C [ (1 + \V^d\^ + \Vd\^ 



+ \u\^\Vu\^ + + |Vdr|VM|" + \u\^\V^d\^ + iVd^dxdr. 

Choosing s = 0, then we obtain 

j (|VMp + iv'dp + \Vd\^){t)dx + j j (p|h|2 + + |V(|Vd|2)|2 

+ |V^d|2+ |Vdt|^) dxdr 


<c + c 



+ I + |Vd|^ + I VdH + \u{^\Vu\'^ + \u\^\V^d\^ 


+ \Vd\^\Vu\^ + \Vdf) dxdr 

t ^ 

C + C [ [{1 + I+ \Vd\^)dxdT + ^ II 2 ^. 


2=1 


Applying Lemma [2.31 to // 2 i(l < i < 5) repeatedly, then it arrives at 


2 r2 

Ihi < sWVWh. + C{e)\\Vd\\^^\\V^d\\l., 


2 r-i 

I 


1122 <e / \\Vu\\jjrdT + C{e) / \\u\\2n\\Vu\\i2dT 


< 


^ I (|lv^i||i. + |||V<i||V 2 <i|||^.)<ir + C(e) / ||«||'!- 7 ’||Vu|||.dr, 


<e / (llx/pi|li. + l|VVr„.)* + C(£) 



(4.21) 


(4.22) 


'0 

2x2 

\\Vd\\^^\\Vmi2 + \\u\\^^\\Vu\\l2]dT, 


1123 <e I ||V^d||^idr + C'(£) n\u\nf\\V^d\\l 2 dT, 


n24<e I \\VuWi,rdT + C{e) I WVdf^rTWVuWl^dr 


2^2 


<■: / (|lx/p*lH. + l|V"<if„.)<ir + C(e) / ||V<i||;ir (II V=d|H= + || V«||i,) dr, 


ih,<E I \\\vdnhdT+c(E) I iivdii'5-iiiv<inii,dT 

Jo Jo 

<E r(||V<i||l. + ||V(|Vdn||i=)<ir + C(T) /‘llVdll^^llVdlll.dr. 


r ' 
1^11 

2||2 
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Substituting // 2 i(l < i < 5) into fl4.22l) and choosing e small enough, then it is easy to deduce 
that 


{\Vu\^ + \V^d\^ + \Vd\^){t)dx + 
+ \V^d\^ + \Vdt\^) dxdr 



(p|h|2 + |Vdnv2d|2 + |V(|Vdp)p 


2t-2 


<C + Cj^ (^1 + \\u\\2.f + II Vd||5^ J (1 + II Vnlli, + II + II VdlliO dr, 

which, applying fl4.ip and Gronwall inequality, completes the proof. 

As a corollary of Lemma [4.21 it is a direct result from (11.Ih More precisely. 

Corollary 4.3. Under the condition fl4.ip . it holds that for 0 < T < T*, 


sup ||(ijL2 + / ||M||^2dr < C. 

0 <t<T Jo 

Now, we give the second important estimate—norm of ||i/pnt||L 2 and ||V^d||x, 2 . 
Lemma 4.4. Under the condition 04.ip . it holds that for 0 < T < T*, 

sup {\\V^t\\l. + \\VM\l^)+ [\\\Vu,\\l, + \\V%\\l,)dT<C. 

0<t<T Jo 

Proof. Step 1; Differentiating (II.ip ^ with respect to t, we get 

putt + pu ■ Vut - div{2p.{p)D{ut)) + VPt 
= —ptUt — ptu ■ Vu — put ■ Vu + dw{2p'ptD{u)) — div(Vdt 0 Vd + Vd Q VdJ. 

Multiplying (I4.25p by Ut and integrating (by parts) over D, it arrives at 


□ 


(4.23) 


(4.24) 


(4.25) 


< C y (plulluillVutl + plnllVunutl + plunv^llutl + plMplVullVuil 

7 

+p\ut\^\Vu\ + |Vd||Vdt||Vni| + |n||Vp||Vn||Vn 4 |) dx = Y.Ihi. 

i=l 

Using (14.7p . Hblder, interpolation, Sobolev, and Young inequalities repeatedly, we get 

Ihl < ||\/p||L°°||u||L6|l\/PUtlU3||VMt||L2 < C'||VM||i2||y^Mt||22||\/pMi|||,6||VnjL2 

< C'll\/puJl2||Vni||22 < C{e)\\^ut\\l2 + e\\Vut\\l2, 

Ih2 < l|p|U-||u||L6||Vn||i3||Mt||L6 < C||VM||i2||VM||L6||VMjL2 
<C(e)||Vn||^,+£||Vn,||i2, 

ddss < ||p|U-||u||L6||V^u||L2||Ui||LS < C'||Vn||L2||V^M||L2||VMi||L2 

<C{e)\\VMh+e\\Vu,\\l2, 
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Ihi < l|p|U-||w|li6||VM||L6||VMji2 < C||VM||^2||VM||j:/l||VMi||L2 

<C{e)\\Vu\\l,+e\\Wut\\l,, 

Ih5 < \\y/p\\L^\\Ut\\L4y/PUt\\L4^u\\L2 < C\\VUt\\L4VPu4l2\\VP'^t\\le 

< C\\^u4l4Vu4l < C{e4^u4l,+s\\Vu4l,, 

Ihe < \\vd\\L4Vd4L4Vu4L2 < c\\vd\\H4^dt\\l4Vd4i4Vu4L2 

< C{e,v)\\Vd4l,+4Vd4lr+e\\Vu4l,. 

Substituting < i < 6) into fl4.26p . then we have 

~ j p\ut\‘^dx + ^ j\Vut\^dx 

< Cm\V~P^t\\l2 + IIVnll^O + C(£, v)\\Vd4l. + r^ll+ £||Vn,||i2 (4.27) 

+ (7 y |M||Vp||VM||VMf|cia:. 

Using 04.101) . 04.71) and interpolation inequality, it arrives at 

Mm + ||P|ki < C{\\Mpu 4 l^ + ll|w||Vn||U2 + \\\Vd\Md\\\L2 

< U(||^Mi||L2 + ||n||i6||VM||i,3 + ||Vd||Loc||V^d||L2) 

< Ci\\^u 4 p 2 + ||Vn||i||Vn||i + \\Vd\\m) 

^ C{\\y/pu4L2 + ||V(i||_H-2 + 1) + -||Vm||_h-i, 

which implies 

11^11^2 + ||P||i^i < U(||i/pMi||L2 + ||V(i||H2 + !)• (4.28) 

Combining 04.27P with 04.28p . it arrives at 

< C(e)(||Vp«,|li. + livnili. + 1) + C{e,r,)\\\/d,\\l, + .,11 Vd,||y + £||V«,||i. (4.29) 

+ C / |M||Vp||VM||VMt|dx. 


To estimate the term J |n| | Vp| | Vm| | Vtitldx. Indeed, by using 04.ip . Holder, Gagliardo-Nirenberg 
and Young inequalities, we obtain 

J |M||Vp||VM||VMf|dx < ||M||L6||Vp||L9||VM||^^||VMt||L2 

< C'||Vn|U2||Vp|U9||Vn||y ||Vn||^||Vni|U 2 

g+6 

< C||V■u||^^'\,^||V■u^||L2. 

By Lemma 12.11 it is easy to deduce that 


'u||v[/2,r + ||P||j^l,r < (7(1 + ||Vnt||i2 + ||V^(i||^2). 


(4.31) 


33 
















J.C.Gao, Q.Tao, Z.A.Yao 


9+6 , 


By Young inequality, Combining fl4.30p with fl4.31l) yields 

j \u\\Vp\\Vu\\Vut\dx < C {1 + ||VMt||L2 + \\V^d\\\2)^\\Vut\\L^ 

<C{e){l + \\V^d\\l2)+e\\Wu4l2. 

Substituting fl4.32p into fl4.29p and choosing £ small enough, we obtain 
f P\ut\'^dx+^ [\Vut\‘^dx 


2 dt 


C 


< C(1 + WVputWh + l|V^d||i 2 + liv^dlli.) + C{r,)\\Vd,\\l 2 + v\\Vd,\\l,. 

Thanks to the compatibility condition and (14.7p . we get 

\\V^t\\l2+ [ \\Vut\\l2dT 

do 

<C{t]) + C [ {l + \\^Ut\\l2 + \\V^d\\l2 + \\VWL^)dT + rj [ \\V%\\l2dT. 


(4.32) 


(4.33) 


Step 2: In order to the control the term || V^(i|||^ 2 dr, we should derive some estimate for 
||V^d||^ 2 - More precisely, multiplying fl4.13p by VAdt and integrating (by parts) over f2, we 
obtain 


J \'VAd\‘^dx + J \Adt\'^dx 
d 


d 


(4.34) 


= — / V(m ■ Vd - IVdl^d) ■ VAd dx - j — [V(m • Vd - |Vd|^d)] ■ VAd dx. 


dt 


To give the second term of the left hand side of fl4.34p hrst. Indeed, it is easy to deduce that 

- J ^ [V(m ■ Vd - I Vdpd)] ■ VAd dx 

<C f (|Vut||Vd| + |Vu||Vdt| + lutllV'dI + lullV'dtl + iVd^llV'dl + iVdllV'd^l ^4_35) 


+ |Vd||V^d||d 4 | + iVdHVdil) |VAd|dx = 5 ^ 1 / 4 *. 

i=l 

By using 04.71) . 04.1UI) . Holder, Sobolev and Young inequalities repeatedly, we obtain 

// 4 I < C'||Vd|Uoo||VUt|U2||VAd|U2 < C'||Vd||^2||V='d|U2||VUt|U2 
<C(d)(||V=^d||i2 + l)||V=^d||i2 + d||Vut||i., 

//42 < C||Vu|U3||Vdi|U6||VAd|U2 < C'||Vu||^,i||V3d|U2||Vdi||^,i 

<C{e)\\Vu\\U\^^d\\h+4'^dt\\l, 

< C{e){\\^u,\\l2 + ||V=^d||i2 + l)||V^d||i2 + e||Vd*||^„ 

//43 < ||ut||L6||V"d|U3||VAd|U2 < ||VuJz.H|V'd||j,i||V=^d|U 2 

<C'(d)(l + ||V^d||i2)||V3d||i2 + d||Vu,|li2, 


34 





















Strong Solution to the Density-dependent Incompressible Nematic Liquid Crystal Flows 


Ih, < C\\u\\Lo.\\V%\\L4VAd\\L2 < 

< m\^put\\l. + ||V=^ci||i. + l)\\VM\l^+e\\W^d4l., 

//45 < C\\Vdt\\L4^^d\\L3\\VAd\\L2 < C\\Vdt\\H4^"d\\H4'^^d\\L2 

<C'(£)(l + ||V=^rf||i.)||V3rf||i. + £||Vrf,||^„ 

//46 < C\\Vd\\L4VmLe\\d4L4VAd\\L2 < CH Vd||j,i || ll v3ci|U2 

< C'(e)(l + \\V^d\\l 4 \\VM\h + e{\\Vd 4 h + 1), 

Ihr < C\\Vd\\L4\V%\\L2\\VAd\\L2 < C\\Vd\\H4^''d\\L2\\V%\\L2 

<C'(e)(l + ||V=^ci||i.)||V=*d||i.+e||V^d*||i., 

//48 < C\\Vd\\le\\Vd4Le\\VAd\\L2 < C|| Vrf|| || Vrf* || || 1^2 

<Cie4Wmi.+e\\Vd4l,. 

Substituting 1 / 4 * (1 < i <8) into (I4.35p . then we obtain 

- / ^ [V(« ■ Vrf - I VrflU)] . VAd dx 

< c(e,s)(i + II vp«,iii= + iivnii|,)iivniii,+£||v<i,fH. + d|iv«,iii=. 


Substituting (I4.36p into (I4.34p and integrating the resulting inequality over (0,t), it arrives at 
\'V^d\‘^dx + [ [ \'V‘^dt\‘^dxdT 


1 

2 


<C + C{e4) / 4 + \\VP^4h + W^Mlhn^WL^dr + e / \\V%\\i.dT 
Jo Jo 

+ dj^ \\Vu4l2dT + C J{\Vu\\Vd\ + \u\\V^d\ + \Vd\\V^d\ + \Vd\4\VAd\dx. 
Choosing e = ^ we obtain 


\V^dfdx+ \V^dt\'"dxdT 



2j |2, 


<C + Ci 6 ) / {i + \\^u 4 i 2 + \\V^d\\i24^M\4dT + 6 \\Vu 4 i 2 dT 
Jo Jo 

+ C j (|Vn||Vd| + \u\\V^d\ + iVdllV^dl + \Vd\4\VAd\dx 

= C + C{6) [\l + \\^u4l2 + IIV=^d||i.)||V=^d||i.dr + 5 r IIVutlli.dr + ^ Ih,. 
Jo Jo 


(4.37) 


By using the Holder, Gagliardo-Nirenberg and Young inequalities repeatedly, we get 


Ihi < C||Vd|Uoc||Vn|U2||VAd|U2 < C||Vd||l.||Vd||^,||Vn|U2||VAd|U2 

< C(1 + liv^diuoiliv^diu^ < C{e) + e\\vmi 2 , 

//52 < C'llnlUellV^dlUsIlVAdlU^ < C||Vn|U2||V2d||il|V2d||i||VAd|U2 

< C(1 + ||V3d|U2)5||v3d|U2 < C{e) + ellV^dlli., 
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Ih3 < C\\Vd\\Le\\V^d\\L3\\VAd\\L^ < C\\Vd\\H4'^^d\\l\\V^d\\l,\\VAd\\L2 

< ^(1 + < C{e) + e\\vmi„ 

IhA < C\\Vd\\l4VAd\\L2 < C\\Vd\\l^\\V^d\\L2 < C{e)+e\\V^d\\l2. 

Substituting 1 / 5 *(1 < i < 4) into fl4.37p and choosing e small enough, it arrives at 

l|V'd||i2+ f \\V^dt\\l,dT 

Jo 

<C + C{5) f {l + \\^ut\\l2 + \\V^\l,)\\VM\hdT + 5 [ WVutWl^dT. 

Jo Jo 

Adding fl4.38p to fl4.33p and choosing S and rj suitably small, we obtain 

iW^/p^tWh + + AllVutlli. + \\V%\\h)dT 

Jo 

<c + c [\i + \\4pu4l, + \\vmi4{i + \\vmi4dT, 

Jo 

which, together with Gronwall inequality and (14.7p . completes the proof. 

Finally, we derive the high order estimate for the strong solution {p,u,P,d). 
Lemma 4.5. Under the condition fl4.ip . it holds that for 0 < T < T*, 
sup iWptWLi + \\u\\h^ + \\P\\h^ + W'VdtWi^ 

0<t<T 

+ [ + IMnllL2 + ||V^(i||^2) dt < C. 

Jo 

Proof. By fl4.28p and (14.311) . it is easy to deduce 

II'*^I|h 2 + ||-P||i/l + [ (||■*^||w2.^ + ll-^llwl-’') — ^5 
Jo 

which, together with (11.ip ^. yields 

llptllL'j < ^ c'||m||j^ 2 ||Vp||l9 < c. 


(4.38) 


□ 


(4.39) 


(4.40) 


(4.41) 


By (I4.13p . (14.251) . (I4.40p and Sobolev inequality, we obtain 

WVdtWi^ = ||VAd +V(|Vdpd-n-Vd)|U2 

< Ci\\V^d\\L2 + |||Vn||Vd||U2 + |||n||V"d||U2 + ||Vd||i 6 + \\\Vd\\VM\\L 4 

< C{\\V^d\\L2 + \\Vu\\L2\\Vd\\H2 + \\u\\H2\\V^d\\L2 + llVdll^i + \\Vd\\H2\\V^d\\L2) 

< c. 


Differentiating (II.ip ^ with respect to t, it arrives at 

dtt-Adt = {\Vd\^d-u-Vd)t. 


(4.43) 
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Taking estimate to fl4.43p . by virtue of fl4.24p and fl4.40p . we obtain 


/ Wduh^dt 

Jo 

<C [ {\\Adt\\L 2 + |||Vdi||Vd||U 2 + \\\\/d\^\dt\\\L^ + \\\Ut\\Vd\\\L. + \\\u\\\/dt\\\L^)dt 
Jo 

<C f {WV^dth^ + \\Vdt\\L4^d\\H^ + \\dt\\L4^d\\jj2 + \\VUt\\L4^d\\H^ 

Jo 

+ ||M||i/2|| V(it||L2) dt 

<C [ i\\V^dt\\L^ + \\Vut\\L2 + l)dt < C. 


(4.44) 


Taking V operator to fll.ip ^. we get 


-V^Ad = V^{\Vd\‘^d-u-Vd-dt). (4.45) 

Applying estimate to d with boundary = 0 and (14.451) we get 

\\\/^d\\L2<C{\\V^Ad\\L2 + \\V^d\\m) 

< C'dlV^dVdpd -u-Vd- dt)\\L 2 + 1) 

< C {\\V%\\\l2 + IllV^dldU^ + \\\Vd\^\V^d\\\L2 + \\\V\\\Vd\\\L2 

+ \\\Vu\\Vm\L2 + \\\u\\V^d\\\L2) (4.46) 

< C {WV^dtlh^ + \\V^d\\l, + \\Vd\\l4V^d\\L^ + \\Vd\\H4^ML^ 

+ II^'*^IIhi II + ||m||// 2|| V^(i||£2) 

^Cdiv^ddiu^ + i), 


which means 

V^dWl^dt < C, 

which completes the proof. □ 



After having the Lemmas 14.1114.51 at hand, it is easy to apply the Theorem II.II to extend the 
strong solution {p,u,P,d) beyond time T*. Therefore, we complete the proof of Theorem 11.21 


5 


Proof of Theorem 


1.5 


In this section, we will also give the proof of Theorem 11.51 by contradiction. Assume 0 < 
T* < oo to be the maximum time for the existence of strong solution (p, m, P, d) to fll.ip - fll.4p . 
Suppose that fll.lOp were false, that is 

- lim (||Vp||loo(o,T;LP + ||V(i||L'>( 0 ,T;L-)) < oo. (5.1) 
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Lemma 5.1. Under the condition fIS.ip . it holds that for 0 < T < T*, 


V^dfdxdr < C, 


(5,2) 


where C denotes generic constants depending only on Q,Mi,T* and the initial data. 
Proof. By virtue of \d\ = 1, we have d ■ Ad = |V(ip. Then, by fl4.2p . we obtain 


f f \Ad\^dxdT = f [\Vd\^dxdT+ f 

'0 J Jo J Jo ■ 

rT r 


I Ac? + \Vd\^d'^ dxdr 


< J J \Vd\^dxT + C. 

By Lemma 12.31 and fl4.2p it is easy to deduce 

j\Vd\*dx^ J |V<ip|V<ip(ii<e||V<i|||,i+C(£)||V(i||jJ||V<i||'i2 

<£|ivniii. + c(£)(i + iivd|i)|). 

Applying the elliptic regularity with Neumann boundary, it arrives at 


rp rp 

f \\V^d\\l,dT<C [ (||Ad||i. + ||d||^0^r 

'0 Jo 

i-T r 

72 J||2 


rp 

<e I \\V^d\\l 2 dT + C{e) I (1 +||Vd||^)dr 

Jo Jo 

<e [ \\vmUT + C{e) [ (1 + llVdlliJdr 
Jo Jo 

fT 

72 J||2 


/ \\\J^d\\i,dT + c, 


'0 


where, r and s satisfy the condition fll.lll) . Choosing e = i, we get 


V^dfdxdr < C, 


(5.3) 


which completes the proof. 


Lemma 5.2. Under the condition fl5.ip . it holds that for 0 < T < T*, 

sup_(||vti||i. + ||vd||l. + ||vMl!») 

+ |V<i|^|Vn|^ + + IViitP) dxdr < oo. 


0<t<T 


□ 


(5.4) 
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Proof. By fl4.2ip . it is easy to deduce, for any 0 < s < t < T, 

j (|Vu|2 + |V"d|2 + \Vd\^){t)dx + j j {p\u? + + iVdVdHl" 

+ |V^dd + |Vdtd) dxdr 

<C + C [ (|Vu|" + iV^dl" + \Vd\^){s)dx + C [ (1 + ||M||ioo + ||Vd||ioo) 


x(l + ||Vu||i. + ||V^d||i. + ||Vd||i4)dr. 


Let 



(p|ud + |Vdnv2d|2 + |V(|Vdd)p 


A{t) =e+ (iVup + iV^dl" + \Vd\^){t)dx + 

J JO 

+ |V^dd + |Vdtd) dxdr, 
then we deduce from fIS.Sp that 


A{t) < CA{s) + C (1 + \\u\\l^ + II Vd||ioc)A(r)dr, 
Jo 

which, together with the Gronwall inequality, gives directly that 


A{t) < CA{s) exp 


C {1 + \\u\\l^ + ||Vd||ioc)dr 


Let 


$(t) = e + sup (IIVulli. + IIV'dili. + IIVd||i4) + / / {p\u\^ + iVhldV^hd + |V(|Vdd)p 
0<T<t Jo J 

+ |V^dd + iVdtd) dxdr, 

then we have ^ 

4>(T) <G4>(s)exp C [ (1 + llwllioo + ||Vd||ico)dr 


(5.6) 


Now we control the term + ||M|lioo + || V(i|||oo)dr. Indeed, by the Lemma 1X41 we have 


C (l + ||M||ioc + ||Vd||ioc)dr 


< c 


(T s) + ^||M||j;^2(g + II Vd|| 2 , 2 (g (in (e + ||m|| 2 , 2 (sJ'. 14 /i, 3 )) ^ ^ 

+ In (e + II V(i||L2(s^'r.vi/i,3)))] . 

Applying the Sobolev inequality and regularity estimate, it arrives at 

hlPwi.3 < < CiW^iiWl. + |||Vd||V^d|||i2), (5.8) 

||Vd||2^i.3 < C'||Vd|p^2.2 < G(1 + ||V2d||i2 + llV^dlli.). (5.9) 
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Substituting fl5.8p - fl5.9p into (15.71) yields 

(1 + \\u\\l^ + l|Vd||ioc)dr < + 

which, together with (15.6p . gives directly that 


Thanks to (14.2p and (15.21) . we can choose s closes enough to T* such that 


lim C 

T^T* 



2 

L^{s,T-,m) 





which, means 


4>(T) < C$(s)2 < cx). 


Thus, we complete the proof. 


□ 


Remark 5.3. Unfortunately, we cannot derive the bound, just depending on the initial data, 
for (15.4p owing to the technigue used here. However, the bound is unform with respect to time 
in (15.4p since s, which closed enough to T*, is fixed in process of the proof for Lemma [Ql 
Thus, we can rewrite (15:41) as 


sup 

te[o,T] 


(|Vn|2 + iV^dl^T |Vd|^) dx 


+ j j {P\M^ + + \Vdt\^ + \V^d\^) dxdt < C{s), 

where and in what follows, C{s) denotes generic constants depending not only on fi, Mi,T* and 
the initial data, but also on the data that is fixed on time s. 

Similar to Lemma 14.31 as a corollary of Lemma (14.2p . we have the following estimate. 

Corollary 5.4. Under the condition (15.ip . it holds that for 0 < T < T*, 

rT 


sup ||dt||L 2 + / \\u\\ff 2 dT < C{s). 
0<t<T Jo 


Now, we derive the high order estimate ||A/p'Ut||L 2 and ||V^(i||i, 2 . 


Lemma 5.5. Under the condition (15.ip . it holds that for 0 < T < T*, 


sup (ll\/P'Ut ||^2 + ||V^d||^ 2 ) + f (||VMt ||^2 + IIV^dt||^ 2 )(ir < C{s). 

0<t<T Jo 


(5.10) 


(5.11) 
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Proof. Recall from fl4.29l) . we have the following estimate 

ITt / 5 

< C(e)(|| + llV’dfi. + 1) + C{e,r,)\\Vdt\\l, + +£||V«,|li. (5-12) 



By fl5.ll) . Holder, Sobolev ang Young inequalities, we obtain 


j luWVpWVuWVutldx < ||M||^^||VM||^^||Vp||L9||VMi||L2 

< C||VM||i,2||VM||Hi||Vp||L9||VMt||L2 
<C{e)\\Vu\\l,+e\\Vut\\l, 

^ C'(||i/p-u||^2 + ||V^(i||^2 + 1) + £||V-Ui||^2. 

Substituting fl5.13p into fl5.12p . we obtain 

l±J Plupdx + y Wupdx 

^ C{£){\\^/put\\L2 + IIV^d||^2 + 1) + C(£,?7)||Vdt||^2 + hllVdtll^i + £||V'Ut|||2. 
Thanks to the compatibility condition, after choosing e small enough, we get 

W'/p’^tWi^ + [ ||V-Ut||L2'^''' 

Jo 

<C{r,) + C [ {l + \\Vput\\h + \\Vmi2)dT + rj [ WV^Wl^dr. 

Jo Jo 

Recall from fl4.38p . we have the estimate 

l|V'd||i2+ r||V'^^tlli2dr 

Jo 

<C + C{5) f (l + ||VP«Ji2 + ||V3d||i2)||V3d||i2dr + <5 / ||VMi||i2dr. 

Jo Jo 

Adding fl5.15p to fl5.16p and choosing e and 6 small enough, we have 

iW^/p^tWl. + ||V=^d||i2) + r(||V«,||i2 + ||V^d*||i2)dr 

Jo 

<C + C [\l + IIVd^tllL + IIV3d||i2)(l + IIV3d||i2)dr, 

Jo 

which, together with the Gronwall inequality, complete the proof of the lemma. 


(5.13) 


(5.14) 


(5.15) 


(5.16) 


□ 
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As a corollary of Lemma 15.51 it is easy to deduce the following estimate 


sup (llptllL'j + 11 ^ 11^2 + ||P||//i + ||V(it||L 2 ) 

0<t<T 


+ [ (ll'^^ll+ IIV‘^(i||^2) dt < C(s). 
do 

Therefore, having all the estimates at hand, it is easy to extend the strong solution beyond 
time T*. Thus, we complete the proof of Theorem 11.51 


6 Proof of Corollary 11.71 

In this section, we will give the proof Corollary 11.71 Indeed, let (p, u, P, d) be the solution of 
fll.lll - fll.4p . we will derive some maximal principle for the direction field. 

Lemma 6.1. For some given constant = 1, 2), then we have the following maximal prin¬ 
ciple: 

i})If 0 < dpj < doi < 1, then 0 < dg, < dj < 1 for any i = 1,2; 

(ii)// —1 < doi < —dgj < 0, then —1 < d* < —dgj < 0 for any i = 1,2. 


Proof. Since (i) has been proved in [3l|, we only give the proof of (ii). Indeed, letting Vi = di+dgj 
and = max{14, 0} > 0, then we obtain 


dtVi - AC, = IVd|2(C - d^f) - u ■ VC. (6.1) 

Multiplying flb.ip by C"*" using the Neumann boundary condition = 0, we get that 

fjfwrrwh + II vr+iii, = j [|v<i|"(i/i - < 4 ,.) - u ■ vv4] • v^dx 

= J\Vd\‘\V*\^dx - j\Vd\%,V*dx + ly dwu\V+\^dx 

< j\^df\vr\^dx < ||V<i|||„||l4+||i,. 

Hence, if we applying Gronwall inequality to the above inequality, we get 

Knmh < IIW(0)||i.exp (^C J‘WVdWl^dr'j =0, 

which implies — 1 < d, < — dg, <0. □ 


Before starting our main result in this section, we recall the elliptic estimate 

||VV||i.<C'i(||A/||i2+ 11/11^0, (6.2) 

and the Gagliardo-Nirenberg inequality 

l|v/||i4<cci|vVlli.||/llioc + ||/||i.), (6.3) 

where the constants Ci and C 2 independent of the function /. 
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Lemma 6.2. For any i{i = 1,2), if the i—th component of initial direction field doi 
the condition 

0 < d^Qi < doi < 1 or — 1 < doi < —< 0 , 
where d^^ is a constant and is defined by 


dioi > 1 


1 

2C1C2' 


where the constant Ci and C 2 are defined in fl6.2p and fl6.3p respectively. Then we 
following estimate 

\\vd\\l^dt < a 

Proof. We only give the proof for the case i = 2. 

Indeed, if 0 < (1^2 < do 2 < 1, then applying the maximal principle, we obtain 



0 < do 2 < <^2 < 1 for any f > 0 and any x G fl. 
By virtne of \d\ = 1, then 


\\d 6211^00 ^ 2(1 ^ 02 ), 

where 62 = (0,1). Snbstituting fl6.4p - fl6.5p into fl 6 . 8 p . we get 

CiC*2||d — e2|l£oo < 1. 


Combining fl6.2p with fl6.3p . we get 

\\V^d\\l.<C,{\\Ad\\l. + \\d\\l,) 

<C,{\\Vd\\l + \\Ad+\VdMl, + \\drH,) 

< CiC2\\d - e2\\U\\V^d\\l, + Cl (C2||d - e2\\U + II+ iVdl^dHi^ + ||d||^i) , 

which, together with fl 6 . 8 p . gives 

II < C {\\d - e2\\U + II Ad + I Vdpdili. + lldll^i) . 

In view of the basic energy inequality fl4.2p . we have 

V^d\\l,dt < C, 


\\Vd\\%dt < C. 

is simple fact if we replaced the function 
we complete the proof. 



which, together with (16.dp . implies 



As for the case — 1 < do* < —dpA < 0 , it 
d + 62 in fl 6 . 8 p - fl 6 .lip respectively. Hence, 


satisfies 

(6.4) 

(6.5) 

have the 

( 6 . 6 ) 

(6.7) 

( 6 . 8 ) 

(6.9) 

( 6 . 10 ) 

( 6 . 11 ) 

( 6 . 12 ) 

' - 62 by 

□ 
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Choosing r = s = 4 in fll.lip . if the maximal existence of time T* < oo, we have 

{\\'^P\\L°°iO,T-,L<i) + ||V(i||L4(0,T;L4)) = CO, 

which, together with fl6.12l) . gives immediately 


^lhnJ|Vp||Loo(o,T;L^) = oo. 

Therefore, we complete the proof of Corollary 11.71 
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